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Abstract 

The weak coupling expansion is applied to the single flavour Schwinger model with Wilson 
fermions on a symmetric toroidal lattice of finite extent. We develop a new analytic method 
which permits the expression of the partition function as a product of pure gauge expectation 
values whose zeroes are the Lee- Yang zeroes of the model. Application of standard finite-size 
scaling techniques to these zeroes recovers previous numerical results for the small and moderate 
lattice sizes to which those studies were restricted. Our techniques, employable for arbitrarily 
large lattices, reveal the absence of accumulation of these zeroes on the real hopping parameter 
axis at constant weak gauge coupling. The consequence of this previously unobserved behaviour 
is the absence of a zero fermion mass phase transition in the Schwinger model with single flavour 
Wilson fermions at constant weak gauge coupling. 
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Introduction In lattice gauge theory, there has been considerable discussion on the phase struc- 
ture of gauge theories with Wilson fermions. A system of free Wilson fermions exhibits a second 
order phase transition at (f3, k) = (oo, l/(2d)), and k being the usual inverse gauge coupling 
squared and hopping parameters respectively and d being the lattice dimensionality. Recent dis- 
cussions concern the extent to which this phase transition extends into the {(5, k) plane, the ex- 
pectation being that there is a line of phase transitions (the chiral limit) extending to the strong 
coupling limit (3 = 0. In the simplest model to which this applies, the lattice Schwinger model 
U, there is also a second order phase transition in the strongly coupled limit Q. This critical 
line, which is expected to connect these two phase transitions, is also expected to recover massless 
physics and therein lies the importance in determining its nature and position. 

While earlier numerical results casted doubt on the above scenario ||, more recent ones support the 
hypothesis that the critical line extends over the whole range of couplings |4], B |6| . To our knowledge, 
only two groups have attempted to determine the phase diagram in the lattice Schwinger model. 
While in rough agreement regarding the location of the critical line, they differ in their conclusions 
regarding its quantitive critical properties. Using Lee- Yang zeroes 0, finite-size scaling techniques 
H as well as PCAC relations, the results of ||] support the free boson scenario, where the model 
lies in the same universality class of the Ising model, with correlation length critical exponent v = 1 
and the specific heat (chiral susceptibility) exponent a = 0. This is not in agreement with Q where 
finite-size scaling of Lee- Yang zeroes on larger lattices provides evidence for v = 2/3, a = 2/3. 

The precise analytical determination of the phase structure in the weakly coupled regime is the 
primary motivation for this work. We contend that the Lee- Yang zeroes in the Schwinger model 
with fixed weak gauge coupling display very unusual size-dependent behaviour, do not accumulate 
on the real hopping parameter axis and may not, in fact, lead to a phase transition. The behaviour 
of these zeroes is, however, such as to mimic the appearance of a phase transition when a finite-size 
analysis is restricted to small or moderate lattice sizes. 

Lattice Schwinger Model We consider a finite d = 2 dimensional lattice with spacing a and 
N sites in each direction. For the fermion fields, we impose antiperiodic boundary conditions in 
the temporal (1-) direction and periodic boundary conditions in the spatial (2-) direction. Lattice 
sites are labelled = n^a where = —N/2, . . . , N/2 — 1. With these mixed boundary conditions 
the momenta for the Fourier transformed fermion fields are k^ = 2irk fl /Na, where k\ = —N/2 + 
1/2, ... , N/2 - 1/2 and k 2 = -N/2, N/2 - 1. 

The action for gauge invariant lattice QED with Wilson fermions is Sqed = Sa[4>] + Sj^^ [0, ip, ip], 
where Sp^^ [<p, iJj, tJj] = J2 m n ~ip{m)M^ w ^ (m, n)tp(n) and 

JlfW (m, n ) = - - J] {( r - 7„)D' /1 (m) Wfl,n + (r + T^W^} • (1) 

This Wilson fermion matrix can be written in terms of free and interacting parts as M^ w ^ = + 
M( mt ) where is given by = 1 in (|l]). Here U^{n) = expi</> /i (n) = exp (ieQaA^n)) is the link 
variable. The Wilson parameter, r, is henceforth set to unity. The partition function is Zqed = 



fVUVt/jVipexp({-S G - 4^}) oc fVUexp{-S G }detM<W oc (detikfW), the integration over 
the Grassmann Fermionic fields having been performed. Here (O) is the pure gauge expectation 



where Up 



value of the quantity O and the pure gauge action is Sg[4>] = PJ2p 1 ~~ \ {Up + Up) 
is the usual product of link variables around a fundamental plaquette and where (3 = l/e^a 2 . 
In the weak coupling expansion we employ the Feynman gauge for the calculation of pure gauge 
expectation values. (There, for large enough N, the gauge propagator in momentum space is 
given approximately by (<j>ij,(p)<j>v(—p)) = 5 P ,v{N d a 2d /2(3) Y^ P , P ^o U ~~ cosp p a) with the infra-red 
zero-momentum mode excluded as standard ||.) 

For P = oo, the only gauge configurations which contribute have S@ = 0. We assume for the follow- 
ing analysis that all gauge configurations satisfying this condition belong to a single equivalence class 
under the group of gauge transformations. Since the fermion determinant is gauge invariant under 
the same group of gauge transformations it is sufficient to evaluate the full partition function on this 
single configuration to determine the behaviour of the system at j3 = oo Therefore the partition 
function is simply proportional to detM^ ). Employing lattice Fourier transforms, the free fermion 
partition function can be written in terms of its dN d eigenvalues by solving the d x d dimensional 
eigenvalue problem of the form M^(p)\Xa <(p)) = ^a\p)\Xa\p)} where a = 1, . . . ,d. In d = 2 



dimensions, the resulting eigenvalues are X a (p) = 1/2k — J2 P =i cosppa + j(-)°yEp=i sin p p a. 

With mixed boundary conditions, the eigenvalues in the free fermion case are either two-fold 
(p2 = or —N/2) or four-fold degenerate. In this case, and with r\ = 1/2k, the Lee-Yang 
zeroes are f]a\p) = Y? P =i r cos Pn a ~ *(~) a \/Ep=i sinp p a. The lowest zero with finite real part 
in k corresponds to p = (±1/2,0) and for a symmetric lattice is r/i°^(±l/2, 0) = 1 + cos(7r/iV) — 
£(— )°sin (tt/N). Pinching of the positive real finite hopping parameter axis occurs at the massless 
point k c = \j2d and application of finite-size scaling [|| to the imaginary parts of the zero gives the 
critical exponent v = 1 (a = then follows from hyper scaling). Therefore the free fermion model 
is in the same universality class as the Ising model in two dimensions and describes free bosons. 
In the presence of a gauge field the position and nature of this phase transition will change and a 
quantitive investigation into the nature of that change is the central aim of this paper. 



The Weak Coupling Expansion The weak coupling expansion is formally obtained by expand- 
ing the link variables U p (n) as a power series in eo- The weak coupling expansion of the matrix 
M( int ) in ([[]) is M( int ) = MW + + O(ejj) where is of order eft. The fermion determi- 

nant detMW = detM(°) x det MC'^MW = det M<W exp tr In (1 + M^^M^) can then be 
expanded giving the following 'traditional' additive form for the weak coupling expansion |jLl| |: 



(detM("0) _ d f* U 1 d f % (2) 

in which U = (M^ t] ) and % = tji = (M^ nt) M^ t] ) - (M^ nt) MjJ nt) ). Here the index % generically 
stands for the combination of Dirac index and momenta (a,p) which label fermionic matrix elements 
so that M-j nt ^ represents (Xa\q)\M^ nt \q,p)\X^\p)}. We note at this point that the expansion is 
analytic in rj with poles at rj = ry| . We note further that tj = (M^ nt ^) at order e is proportional 



to the pure gauge expectation value of the gauge field 4>^ and is zero. 

An alternative formulation of the partition function may be obtained by writing the Wilson fermion 
matrix as M( w > = r/ + H where H is the hopping matrix. The fermion determinant det M^ W ' = 
det(r/ + H), is a polynomial in r\ since for finite lattice size these matrices are of finite dimension. 
Therefore the pure gauge expectation value of the fermion determinant is also a polynomial in r] and 
as such may be written (det M^ w ^) = Y\f=i (rj — r]i). Here rfc represent r] a (p) and are the Lee- Yang 
zeroes in the complex 1/2k plane. We may write a 'multiplicative' weak coupling expansion as 

detMW tU-.f) Ml V-VP) ' ^ 

where Aj = rji — rjf 1 ^ are the shifts that occur in the zeroes when the gauge field is turned on and 
are the quantities to be determined. Note again that (||) is analytic in rj with poles at r]\ . Since 
the two expressions @ and (g) must be equal, the residues of the poles must be identical. 

Let ti = tf ] + ...,Uj = t\f + . . . and Aj = r/f ) + r)f ) + . . . where r/f ) is 0(e) and tf\ t\f and 
rff^ are 0(e 2 ). Assume the n* h free eigenvalues are degenerate for j = 1, . . . , D n where D n is 2 
or 4 and An/ = A n . At this stage we make no assumptions regarding the other eigenvalues. The 
traditional weak coupling expansion (Q) and its 'multiplicative' counterpart (||) have single poles 
at 7] = r]n \ the residues of which are easily calculated to 0{e 2 ). Equating the two residues order 
by order in e gives 



D n D n D n 

j=l j=\ j=i 



-4? + E — ^ 



_(0) _(0) 



(4) 



In this latter equation, An is the average of two zeroes which in the free fermion case were two-fold 
degenerate and is 0(e 2 ). Equating the residues of the corresponding double poles gives no extra 
information at this order. 

In the non-degenerate case D n = 1, (||) is sufficient and necessary to fully determine the shifts in 
the positions of the zeroes r)n) and rjnj . This then amounts to a new method to uniquely express 
the partition function and all derivable thermodynamic functions in multiplicative (as opposed to 
the traditional 'additive') form. Thus, we have expressed the pure gauge expectation value of a 
product of fermion eigenvalues as a product of pure gauge expectation values. 

In the -D n -fold degenerate case, @ gives for the expansion of the 'multiplicative' weak coupling 
expression (g) to 0(e 2 ), 

( A, \ t (2) tf) AjA • 

ni 1 (oj ) = 1 + £7(07 + £ : (o) ; : (o) 7W\ + ^ TWTM + ' ' ' ' ^ 



This is identical to the traditional weak coupling expansion (||) to 0(e 2 ) if 



Setting Tj = T)n +e where e / and separating out the contributions coming from the D n -fold 
degenerate cases rij now gives that 

D n D n D„ 

E(^) 2 = E E *g^+ow • (7) 

Now the 'multiplicative' and 'additive' expressions for the partition function coincide up to 0(e 2 ) 
everywhere in the complex 1/2k plane arbitrarily close to any pole and have the same poles and the 
same residues (up to 0(e 2 )) at those poles. The partition function zeroes are rj ni = rjn^ +r]n} + 
where in the free field case the nj th zero is in the re th degeneracy class. In terms of the Dirac and 
momenta labels, we can write the erstwhile 2-fold degenerate lowest zeroes as 

7 la (± Pl ,0)=r ] £\\p 1 \,0)±r ] £\\p 1 \,0) + A£\\p 1 \,0) . (8) 

In this way, the shifts in the positions of the erstwhile two-fold degenerate zeroes are also determined 
to C(e 2 ) and the 0(e 2 )-shift is the shift in the average position of the two zeroes while their relative 
separation is 0(e). The calculation of Lee- Yang zeroes is now reduced to the calculation of the 
pure gauge expectation values A a (p) which we do in Feynman gauge. 



Results and Conclusions The main plot in Figure 1 is a finite-size scaling plot for the lowest 
zero for (3 = 10 and N = 8,10, ...,62 coming from our weak coupling analysis (circles). The 
corresponding data from the numerical analysis of Q are also included (crosses) . The two lines are 
the linear fit to the lowest zeroes for lattice sizes 16, 20, 24 of || and a fit to the second zero. These 
yield slopes —1.5 and —1.4 respectively, corresponding to v ~ 2/3. The insert contains a similar 
plot at (3 = 5 where the circles are our results for N = 2, 4, 42 and the squares are the results of 
0] for the lowest zero for iV = 2, 4, 8. The line is the best linear fit to this latter data and yields v 
compatable with 1. Fitting to our small lattice data also gives v ~ 1. 

Therefore, confining the finite-size scaling analysis to the range TV = 2 - 8 yields a slope compatable 
with v = 1 in agreement with 0|. A corresponding analysis with N = 8-24 gives a steeper slope, 
compatable with v = 2/3 in agreement with Q. It is clear from the figure, however, that the 
curve does not in fact settle to a finite-size scaling line. Instead as ./V increases, the lowest zero 
crosses the real 1/2k axis. The second zero exhibits the same behaviour as demonstrated by the 
upper curve in Figure 1. The first two zeroes therefore fail to accumulate and do not contribute to 
critical behaviour 12, 13]. One interprets these zeroes as isolated singularities of thermodynamic 
functions, having measure zero in the thermodynamic limit. 

It would be extremely difficult to detect the isolated nature of these Lee- Yang zeroes from numerical 
data alone. A standard numerical technique for determination of the phase diagram - assuming 
accumulation of zeroes - is to approximate the infinite volume critical point by the real part of 
the lowest zero (then a pseudocritical point) for some reasonably large lattice. Plotted against (5 
this appproximates the phase diagram. In Figure 2 we present such a plot for N = 24 (circles) 
to compare with the results of || (crosses) also at N = 24. The phase diagram of || for the 
Schwinger model with two fermion flavours coming from a separate PCAC based analysis is also 
included (squares) for comparison. 
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Figure 1: The finite-size scaling behaviour of the imaginary parts of the first two Lee- Yang zeroes 
for (3 = 10. The symbols o are our weak coupling expansion results. The numerical results of || 
for N = 8, . . . , 24 are x and the lines are best fits to these points. The insert is a similar plot for 
P = 5. There the □ are the lowest zeroes of Q. 



Now, the 0(e) shift in the free field lowest zero is, in fact, ^(±1/2, 0) = ±(i/(Na) d ) y J{\^ (1 .(l, 0)|) = 
±i/(2y/~/3N\ sinpx^D) the leading finite-size scaling behaviour of which is AMndependent. The 
crossing of the real axis by the first two zeroes is therefore due to the N dependency of their C(e 2 ) 
average. In the four fold degenerate case, the leading N behaviour of the zeroes is again determined 
by that of their average. The four fold degenerate C(e 2 ) averages of small momentum zeroes also 
cross the real axis and do not accumulate. 

Although we cannot logically exclude the possibility that some zero which was far from the real axis 
in the free fermion case (or some combination of such zeroes) conspires to accumulate on the real 
1/2/c axis in the usual way, such a situation seems unlikely. Nor can we exclude the possibility that 
just as the subdominant terms in the weak coupling expansion (the 0(e 2 ) averages) correspond 
to the dominant N behaviour, so too could higher order terms in the weak coupling expansion 
correspond to superdominant A" effects. This seems again unlikely as the average zeroes coming 
from our weak coupling expansion agree very well with those of (see |]l0| for more details). 
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Figure 2: The would-be phase diagram coming from the real parts of the lowest zeroes at various (5 
at fixed N. The symbols o are our weak coupling expansion results for N = 24 and x are numerical 
results of || also for N = 24. The symbols □ are the numerical results of || for N = 16 and two 
fermion flavours. 



Moreover, superdominant N effects from subdominant weak coupling terms would spell disaster 
for lattice (and indeed continuum) perturbation theory and is an unlikely scenario in this - a 
superrenormalisable - model. 

Although the possibility of existence of isolated singularities and non-accumulation of partition 
function zeroes has been known for a long time jj], 12, 113], this is to our knowledge the first 
instance where such behaviour has been observed. In conclusion, in the free case, there is a phase 
transition precipitated by the accumulation of Lee- Yang zeroes on the real hopping parameter 
axis. In the weakly coupled regime at fixed j3, this accumulation no longer occurs. Instead, the 
movement of zeroes for small and moderately sized lattices mimics phase transition like behaviour. 
As the lattice size becomes large, however, these zeroes move across the real axis, and do not give 
rise to a phase transition. 
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